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ABSTRACT 

We complement a recent ApJ Letter by Luo et al. by comparing the fraction of starburst galaxies 
which are interacting with the overall fraction of interacting galaxies in the nearby galaxy population 
(within 40 Mpc). We confirm that in starburst galaxies the fraction of interacting galaxies is enhanced, 
by a factor of around 2, but crucially we do so by studying a sample of almost 1500 of the nearest 
galaxies, including many dwarfs and irregulars. We discuss how adjusting the starburst definition 
influences the hnal result and conclude that our result is stable. We find significantly lower fractions 
of interacting galaxies than Luo et al. did from their larger but more distant sample of galaxies, 
and argue that the difference is most likely due to various biases in the sample selection, with a 
representative sample of the nearest galaxies, such as the one used here, being the best possible 
representation of a general picture. Our overall conclusion is that interactions can and do increase the 
number of starburst galaxies, and that this result is extremely robust. By far most starburst galaxies, 
however, show no evidence of a present interaction. 

Subject headings: Galaxies: interactions — Galaxies: spiral — Galaxies: star formation 


1. INTRODUCTION 

In current cosmological galaxy evolution models, inter¬ 
actions between galaxies play a particularly prominent 
role in shaping both the dark halos and the baryonic 
structure of modern-day galaxies. Almost equally im¬ 
portant is enhanced star formation in the form of what 
is commonly called starbursts, which not only trans¬ 
forms large amounts of gas into stars in short timescales, 
but also produces substantial changes in the interstellar 
and intergalactic media, through winds and violent late 
stages of stellar evolution. 

Starbursts are related to interactions in the sense that 
the latter can stimulate the former, and many starbursts 
are observed to occur in interacting or merging galaxies. 
There is abundant anecdotal and statistical evidence to 
back up these general statements (see, e.g., Toomre & 
Toomre 1972, Larson & Tinsley 1978, or the review by 
Schweizer 2005), and in particular in the case of the most 
extreme starbursts, such as the luminous and ultra lu¬ 
minous infrared galaxies (LIRGs and ULIRGs), it has 
long been clear that there is a strong statistical connec¬ 
tion with galaxy interactions and mergers (e.g., Joseph & 
Wright 1985, review by Sanders & Mirabel 1996). Such 
extreme starbursts are very rare though, especially at 
the current cosmological epoch—there are no ULIRGs 
within some 70 Mpc of us, for instance. 

It is therefore vitally important to separate the anecdo¬ 
tal from the statistical when referring to the interconnec¬ 
tions between starbursts and interactions. A large body 
of observational and numerical work has led to the gen¬ 
eral conclusions that indeed there is a causal connection, 
but also that statistically the increase in star formation 
rate (SFR) as a result of a galaxy-galaxy interaction or 
merger is limited (e.g., Larson & Tinsley 1978; Bergvall 
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et al. 2003; Kapferer et al. 2005; Di Matteo et al. 2007, 
2008; Ellison et al. 2008, 2010, 2013; Robaina et al. 
2009; Knapen & James 2009; Rodrighiero et al. 2011; 
Saintonge et al. 2012; Moreno et al. 2015). In addition, 
it is clear that most interacting galaxies at present do 
not have an enhanced SFR at all (e.g., Knapen, Gister- 
nas & Querejeta 2015, hereinafter Paper II). This may 
well be due to timescales, in the sense that such ‘qui¬ 
escent’ interacting galaxies may have had an enhanced 
SFR in the past or will have it in the future, but it is 
important to establish observationally, and confirm by 
simulations, the duration and intensity of a typical and 
possible ‘starburst’ phase. 

In this Letter, we report on the results from a study of 
the statistical connections between starbursts and inter¬ 
actions in a representative sample of 1500 of the most lo¬ 
cal galaxies, focusing on the fraction of starburst galaxies 
that is interacting, compared to the fraction of interact¬ 
ing galaxies in the overall galaxy population. This work 
forms an extension, and confirmation, of similar stud¬ 
ies performed on more distant galaxies, often selected 
using the Sloan Digital Sky Survey (SDSS; e.g., Luo 
et al. 2014). It is a vital step towards understanding 
to what extent general conclusions on relations between 
starbursts and interactions, often reached on the basis 
of or influenced by studies of rare but striking objects, 
are applicable to the general galaxy population, consist¬ 
ing of significant numbers of dwarfs and other low-mass 
galaxies. 

2. OBSERVATIONAL DATA, SAMPLE, AND 
METHODOLOGY 

We use a sample of 1478 nearby galaxies {D < 40 Mpc) 
selected from the Spitzer Survey of Stellar Structure in 
Galaxies (S'^G; Sheth et al. 2010). We use a number 
of key parameters, namely the SFR (from a combina¬ 
tion of IRAS 60 and 100 pm fluxes; from Querejeta et al. 
2015) and the SSFR (SFR divided by stellar mass, the 
latter from the dust emission-corrected Spitzer 3.6 pm 
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images from the S^G; also from Querejeta et al. 2015), 
and whether the galaxies are interacting (classes A— 
merging, B—highly distorted, and C—with minor dis¬ 
tortions; from Knapen et al. 2014, hereinafter Paper I). 

As explained in more detail in Paper II, for each sam¬ 
ple galaxy we calculate the enhancement in its SFR and 
SSFR, i?(SFR) and A(SSFR), by dividing these param¬ 
eters by the median values for a control sample. The 
control sample is created for each galaxy individually, 
and consists of all those galaxies which are not interact¬ 
ing (and which also do not have a close companion—a 
further category defined in Paper II but which we do not 
use in the current papeiFj, but which are close in mor¬ 
phological type (within ^ 1 numerical class) and stellar 
mass (within ±0.2 in log^M/Mo)). 

Having calculated i?(SFR) and ^(SSFR), we define 
starburst galaxies as those which have values of E 
above a certain number, in particular A(SFR) > 5 and 
A(SSFR) > 4. We show in the next Sections that these 
choices are reasonable, but also investigate the effects of 
using higher or lower limits, and of limiting the sample 
to certain ranges of stellar mass. 

To study the connection between starbursts and merg¬ 
ers we follow the same approach as Luo et al. (2014), 
namely determining how many of our galaxies defined as 
starburst are interacting (so in our categories A, B, or 
C), and then how many of all sample galaxies are inter¬ 
acting. The ratio between these two numbers, H, the 
starburst boost due to the interactions, is the main pa¬ 
rameter reported in this Letter. For our full sample, we 
find that 18% of our starburst and 9% of all galaxies are 
interacting, leading to a boost factor B — 2.0. We make 
further tests by reducing the sample of galaxies studied 
by imposing that they have a certain minimum stellar 
mass, as described below, but they do not alter our con¬ 
clusions. 

3. RESULTS 

The parameter we are mainly interested in here is 
the enhancement of the fraction of interacting galax¬ 
ies among the starburst population, as compared to the 
overall population of galaxies. For this, it is vital to use 
a reliable and meaningful definition of what a starburst 
is, and in Paper II we argued that a reasonable definition 
is an enhancement of the SFR of a galaxy by a factor of 
5 or more compared to its control population (a factor 
of 4 in SSFR). Using this value for U(SFR)= 5, 18% of 
starburst galaxies and 9% of all galaxies are interacting. 
This implies a boost factor B = 2.0 which is the main 
result reported here (for i5(SSFR)= 4, the numbers are 
22%, 9%, and B = 2.4). 

To investigate how robust these results are with chang¬ 
ing starburst definitions, we show in Fig. how B varies 
with different values of our starburst criteria: the en¬ 
hancement cutoffs E for SFR (left panel) and SSFR 
(right panel). As might be expected, we see that the 
difference in interaction fraction between starbursts and 
the complete sample {B) increases as the starbursts be¬ 
come more extreme (higher E). But what Fig. also 

^ The 138 galaxies in our sample that we classified as having a 
close companion in Paper I (but which are not interacting) have 
the same behavior in terms of SFR and SSFR as the overall sam¬ 
ple, which is why we limited our analysis to the galaxies in our 
interaction classes A, B and C. 


shows very clearly is that the boost factor B is robust, 
and in fact does not vary by more than some 25% (±0.4 
and ±0.6 in SFR and SSFR, respectively) over the whole 
range of starburst definitions probed. And this range is 
rather large: E = 2 hardly discriminates starbursts at all 
from star-forming galaxies, whereas E = 8 selects only 
the most extreme starbursts (for comparison. Table 1 in 
Paper II lists the only 18 galaxies in our sample of 1478 
with E > 10). The conclusion from this test is that our 
result is robust, and does not depend significantly on the 
exact value used to define a starburst, nor on whether 
we use SFR or SSFR for the starburst definition^ 

4. DISCUSSION 

4.1. Starbursts are preferentially associated with 
interactions—Yet most starbursts don’t interact 

The main conclusion we reach is that the fraction of 
starburst galaxies that are interacting is significantly 
higher than in the whole population. The difference is 
a factor of 2, and this enhancement is stable even when 
changing the starburst definition (see Fig.[^. This num¬ 
ber is very similar to that reported in the recent Letter 
by Luo et al. (2014), even though the properties of their 
input sample and the fractions of starburst and interac¬ 
tion galaxies among their sample are very different (see 
next Subsection). 

Our study of a sizeable sample of 1478 of the most local 
galaxies of all types thus confirms that indeed the frac¬ 
tion of galaxies associated with galaxy-galaxy interac¬ 
tions and mergers is significantly higher among starburst 
galaxies than among the general population of galaxies. 
And this in turn is further evidence that indeed mergers 
and interactions are among the main causes for the oc¬ 
currence of starbursts in galaxies. We intentionally write 
here “among the main causes”, as in spite of the fac¬ 
tor of two difference in interaction fraction between the 
starburst and overall samples, the fraction of starburst 
galaxies, 22% when using a reasonable definition of what 
constitutes a starburst, is still modest. The enhance¬ 
ment in this fraction related to interactions and mergers 
is undeniable, but most starburst galaxies, by far, still 
occur in galaxies without any evidence for interactions. 
This may well be related to timescales, with both the 
starburst phase and the visible stage of morphological 
evidence for interactions being relatively short-lived, but 
does mean that caution must be applied when discussing 
starbursts, interactions, and their interrelations. 

4.2. Starburst fractions depend on sample biases 

Although our “boost factor”, of H = 2, by which star- 
bursts are preferentially interacting compared to the con¬ 
trol population is remarkably similar to that found by 
Luo et al. (2014), one very obvious difference between 
their results and ours is that their fractions of interact¬ 
ing galaxies are significantly higher. Luo et al. report 
that around 50% of their starburst galaxies show evident 

We note here, as we do in Knapen & James 2009 and in Pa¬ 
per I in more detail, that other definitions of starbursts can lead 
to significantly different populations of galaxies. We use here the 
definition of SFR and SSFR enhancement, as do Luo et al. 2014, 
rather than alternatives such as, for instance, gas depletion time, 
absolute values of SFR or SSFR, or such absolute values normalized 
by area. 
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SFR Enhancement Cutoff 


SSFR Enhancement Cutoff 


Figure 1. Boost factor B by which the fraction of interacting galaxies is enhanced among starburst galaxies as compared to the non- 
starburst control population, as a function of the enhancement cutoff used to define whether a galaxy is called a starburst: the factor E by 
which the SFR and SSFR of a starburst is higher than those of their control sample. Small labels near each data point indicate the fraction 
of starburst galaxies that occur in interacting galaxies (top) and the fraction of all sub-sample galaxies that are interacting (bottom), the 
plotted data points indicate the ratio between these two numbers. The main conclusions of this Letter are based on reasonable values of 
E(SFR) = 5 and E(SSFR) = 4, which yield B = 2.0 and 2.4, respectively. 
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Figure 2. As Fig. ^ but now as a function of the stellar mass cutoff (lower limit) used to define the input sample. The B values in Fig.^ 
are reproduced as trie leftmost points. All for i7(SFR) = 5 (left panel) and i7(SSFR) = 4 (right panel). 


merger features, against just below 20% for their con¬ 
trol sample. In contrast, we find values of 22% and 9%, 
respectively. 

We suspect that this difference is due primarily to bi¬ 
ases introduced by the sample selection. There are sev¬ 
eral important differences in this respect, including 

• Luo et al. (2014) select galaxies at redshifts from 
0.01 to 0.20, with an apparent magnitude cutoff of 
r = 17.72, whereas our galaxies are within D = 
40Mpc and with uib < 15.5. This translates into 
limits in absolute magnitude of roughly —17.5 for 
our sample, and a range from —17.8 to —23.2 for 
Luo et al. 

• The starburst definition is similar in both papers, 
with SFRs that are ~ 5 times larger than the me¬ 
dian SFR of a control sample. But in practice, 
given the “main sequence” (Fig. 2 in Paper II, for 
instance, or Fig. 1 in Luo et al. 2014) the SFR 
increases with stellar mass, and thus with absolute 
magnitude, so that starbursts will have higher ab¬ 
solute SFRs with higher masses. The starbursts. 


as well as the control galaxies, as defined by Luo et 
al. will thus have much higher SFRs. 

• In both cases, the presence of interactions is partly 
based on visual inspection of optical images. This, 
however, will give substantially different results in 
our case of very nearby galaxies which are well re¬ 
solved in the imaging used, and the much more 
distant galaxies in the sample of Luo et al. (2014). 

Many of these biases are hard to quantify, and it is 
almost impossible to correct the interaction fractions for 
them. But one test that we can easily perform is to simu¬ 
late a different range in stellar mass (and thus in absolute 
magnitude and, in an indirect way, distance) in our sam¬ 
ple. We show the results of this in Fig.[^ which indicates 
how both the boost factor B by whidi the interaction 
fraction increases in starburst galaxies (data points) and 
the interaction fractions among the starburst and control 
samples (small numbers near data points) change sub¬ 
stantially when changing the mass range allowed in the 
sample. We see that the interaction fractions increase at 
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Figure 3. As Fig.|^ but now using an upper limit for the stellar 
mass as cutoff. This highlights the behaviour of low-mass galax¬ 
ies (rather than the high-mass galaxies in Fig. [^. Only the plot 
for i?(SSFR) = 4 is shown because for the low-mass (e.g., dwarf) 
galaxies the SFR is a considerably less useful parameter than the 
SSFR. No significantly different behaviour among the lower-mass 
galaxies is seen. 


higher-mass samples, but, in apparent contradiction to 
the constancy of the boost factor of around 2 between 
the work of Luo et al. (2014) and ours, B drops in the 
highest mass cutoff bin. Here we must caution though 
for small subsample sizes, as our highest-mass subsample 
(M* > 10.5 Mq) contains only 292 galaxies (of which 22 
are starbursts and 270 non-starbursts, 4 and 48 of which 
are interacting, respectively). 

The main conclusion here is, in any case, that in spite 
of these differences in sample selection and methodol¬ 
ogy, and possibly resulting biases, the overall result is 
extremely robust: interactions boost the fraction of star- 
bursts, and do this by a factor of around 2 — 2.5. The 
fact that we reproduce this result from the largest sample 
of the most nearby galaxies now available is very signif¬ 
icant, and lends important further support to the con¬ 
clusion that indeed interactions can and do increase the 
number of starburst galaxies (although most starbursts 
are not presently interacting). 

4.3. Starbursts and interactions in low-mass galaxies 

Fig. [3] is like Fig. but using the stellar mass as an 
upper limit rather than the lower limit of Fig. There, 
we wanted to simulate the lack of low-mass galaxies in¬ 
cluded in samples at higher redshifts, but because our 
sample contains also low-mass galaxies, we can use a sim¬ 
ilar analysis to investigate whether the behavior of those 
galaxies with the lowest masses is statistically different 
from the main population. 

Fig.l^shows that this is not the case: there is no signif¬ 
icant different between the boost factor (by which star- 
bursts are preferentially interacting compared to the con¬ 
trol population) for low-mass galaxies, where the lowest 
stellar masses we can study here are < 10®'^ M©. Of 
course the uncertainties are large because only a few of 
these galaxies in our sample are defined as starbursts 
(only 5 of 367 of the < 10® ® Mq galaxies, 1 and 25 of 
which are interacting, respectively), but the lack of ev¬ 
idence or a different behavior qualitatively agrees with 


what has recently been reported by, e.g., belli et al. 
(2014), Moreno et al. (2015) or Stierwalt et al. (2015). 

5. CONCLUSIONS 

We present a detailed study quantifying to which ex¬ 
tent starburst galaxies are preferentially interacting, us¬ 
ing a sample of almost 1500 of the nearest galaxies for 
which we have reliable information on their SFRs, stellar 
masses, and on whether they are interacting with neigh¬ 
boring galaxies. A crucial difference between this study 
and others in the literature is that we base our results on 
a representative sample of local galaxies, including many 
of relatively low stellar mass. 

We confirm that among starburst galaxies the fraction 
of interacting galaxies is enhanced, by a factor of around 
2, compared to the general population of local galaxies. 
Adjusting the starburst definition—the exact threshold 
by which the SFR or SSFR of a starburst galaxy must 
be enhanced compared to its control sample—allows us 
to conclude that our final result is stable. We find signif¬ 
icantly lower fractions of interacting galaxies than Luo 
et al. (2014) did from a larger but more distant sam¬ 
ple of galaxies selected from the SDSS survey, and argue 
that the difference is most likely due to various biases in 
the sample selection, with a representative sample of the 
nearest galaxies, such as the one used here, being the best 
possible representation of a general picture. Our overall 
conclusion is that interactions can and do increase the 
number of starburst galaxies, and that this result is ex¬ 
tremely robust. By far most starburst galaxies, however, 
show no evidence of a present interaction. 
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